Jet-Evolution in the Quark-Gluon Plasma from 
RHIC to the LHC 

^ ; S. Domdey a , B.Z. Kopeliovich a ' & , H.J. Pirner a 

(j . a) Institute for Theoretical Physics, University of Heidelberg, Germany 

. b) Departamento de Fisica, Universidad Tecnica Federico Santa Maria; 

Instituto de Estudios Avanzados en Ciencias e Ingenieria; 
Centro Cientifico-Tecnologico de Valparaiso, Casilla 110-V, Valparaiso, Chile 

Qh 

Ph| The observed suppression of high-p± hadrons allows different explana- 

tions. We discuss two possible scenarios: In scenario 1, parton energy 
loss from scattering in the hot medium is complemented by final state in- 
^- ■ teractions in the resonance matter. Scenario 2 has an enhanced transport 

parameter q which is fitted to RHIC data. For LHC, the two scenarios 



m 

o 
d 

o 



lead to very different predictions for the nuclear modification factor of 
hadrons. In addition, jet reconstruction allows more specific tests of the 
mechanisms responsible for jet quenching. We calculate the distribution 
of partons inside a jet and find different results for the two scenarios. 



1 Introduction 



Investigations of the hot medium produced in heavy-ion collisions largely suffer from both 
its short life time and the small production rate of hard probes. Proceeding from exper- 
iments at RHIC to LHC, these features will be significantly improved. In particular, the 
increase in energy at LHC will lead to a higher temperature compared to RHIC which cor- 
responds to a longer life time. Secondly, full jet reconstruction provides more differential 
information on the mechanisms of final state interactions compared to spectra of leading 
hadrons. 

The large suppression of leading hadrons as observed by PHENIX and STAR coll. [1] 
(see also j2]) has triggered large interest in final state interactions in the hot medium 
produced in heavy-ion collisions. Theoretical models have focussed mostly on the energy 
loss of fast partons during their propagation through the dense medium. Many discussions 
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1 Introduction 



of the range of mechanisms at work at RHIC and possibly at LHC can be found in the 
literature (see for example [3]). 

In the picture of partonic energy loss, hadronization takes place outside of the hot zone 
and is not affected by the medium. Special features may develop due to the nonabelian 
nature of QCD. Coherent radiative energy loss may be important due to the radiated 
gluon interacting with the medium. This explanation corresponds to a less pronounced 
suppression of heavy charmed quark/hadrons in contrast to experimental results [4]. 

Gluon radiation undoubtedly plays an important role for parton evolution at RHIC 
and the LHC. Typical virtualities of partons are around Q = 20 GeV and Q = 100 
GeV, respectively. The resulting parton showers extend over several fm and evolve in 
the medium. Therefore, a separation of this shower from the propagation of partons in 
the medium seems to be unrealistic. From our perspective, however, there are no strong 
indications for coherence yet, because of the finite size of the medium which blurs the 
characteristic dependence of the energy loss on the square of the path length [5J [7J |9] . 

In our view, it is pragmatic to explore parton evolution interleaving radiation and 
interaction with the medium incoherently. Once the parton dynamics has been settled, 
measurements of the high-p^ hadrons and/or jets can be used as a tool to investigate the 
properties of the quark gluon plasma. The further hypothesis of a strongly interacting 
quark gluon plasma (sQGP) [10] can be tested. This idea underlies models with large 
energy loss parameters q, but has not received so much backing from lattice simulations 
where the quark gluon coupling constant was monitored in the region 0.5 < a s < 1 at most 

mi. 

A typical shower at RHIC evolves from a high virtuality around Q = 20 GeV to a 
low perturbative scale Qq = 1.5 GeV and lasts about r cvo = 2 fm. In comparison, for 
longitudinal expansion with an initial temperature of To = 300 MeV at tq = 0.5 fm, one 
may estimate a plasma lifetime of r c = 3.3 fm [For details on these estimates, see Eq. ( fl8l) 
and Eq. ([19"]) .]. From these two time estimates, we find 

7"evo ^ T c . (1) 

The nonperturbative part of hadronization involves the decay of the prehadrons or reso- 
nances at the preconfinement scale Qq where e.g. the vector meson resonances decay into 
4 — 5 pions. Therefore, it is very probable at RHIC that at the end of the evolution reso- 
nances interact with the remainder of the quark-gluon plasma, i.e. with hadronic resonance 
matter [26]. For T < T c , the equation of state of hadronic matter can well be described 
by a hadron resonance gas with an initial density n res ~ T c 3 [12]. For our purposes, we 
neglect the slow cool-down of resonance matter below T c . Such an interaction may be 
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described by a hadronic theory with cross sections slightly larger than hadronic cross sec- 
tions in vacuum. In general, hadronic resonances have radii which are larger than hadronic 
ground states. Therefore, one expects that gg-resonances have larger cross section than 
the gg-ground states. Because of these large cross sections, newly formed prehadrons are 
absorbed. We think that final state effects of these resonances with the hot resonance gas 
play a decisive role in the observed suppression of hadrons in RHIC experiments. 
We therefore would like to advocate a comparison of two scenarios. 

• Scenario 1 uses an estimate of q = 0.5 GeV 2 /fm for the jet transport parameter at 
RHIC. This scenario underestimates the suppression when only parton energy loss 
is taken into account. Additionally, we consider the absorption of the prehadrons in 
the resonance gas, leading to a similar suppression as observed at RHIC. 

• Scenario 2 considers large parton energy loss solely and tunes up this value to q = 4 
GeV 2 /fm in order to accomodate RHIC data. 

Phenomenologically, we propose to test both scenarios on LHC data. For this purpose, we 
extrapolate both scenarios to Pb+Pb collisions at LHC with y/s = 5.5 TeV, a virtuality of 
Q = 100 GeV and an initial temperature of T = 500 MeV. In the LHC case, the average 
time for the evolution from the initial high virtuality to a hadronic scale Qo is much longer 
than the plasma life time and the size of the hot interaction zone, 

7~cvo > T c . (2) 

Therefore we expect that resonances do not play an important role there. 

Another important feature of LHC experiments is the possibility of full jet reconstruc- 
tion up to very small momenta. Therefore we will also make predictions for jet compositions 
down to very small momentum fractions which govern the jet multiplicity. In both sce- 
narios, the respective temperatures and virtualities control the parametric dependence of 
the input data to the modified DGLAP equations we use. The main uncertainty con- 
cerns the question whether final state absorption is important. Should we describe RHIC 
physics in the pure energy loss scenario 2 or in the mixed scenario 1, where nonperturba- 
tive features of resonance interactions occur together with perturbative parton interactions 
during hadronization? We claim that LHC will tell which scenario is preferred with its 
forthcoming experimental results at the end of this year. 
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2 Formalism and Modified Evolution Equations 



2 Formalism and Modified Evolution Equations 

Hard processes in a heavy-ion collision lead to the production of partons with high virtuality 
and high energy. These partons will successively radiate gluons to reduce their virtuality 
and become on mass-shell. This leads to a parton shower and scaling violations in the 
vacuum jet fragmentation functions D^(x,Q 2 ) described by the DGLAP equations [13] 

9 no, a s (Q 2 ) f X dz 



In Ref. [33], some of the authors have developed a model for the parton cascade that 
includes scattering off the partons in the QGP. In this model, we modified the DGLAP 
equations by an additional scattering term. 

In this section, we repeat the formalism for our further considerations. We consider 
both quarks and gluons as hadronizing partons. Indices on the fragmentation functions 
and the splitting functions indicate the respective parton species. 

The formation of a parton shower does not happen instantaneously, but requires a 
certain time. The relevant time is the time a virtual state needs to evolve in virtuality 
from Q 2 to Q 2 + dQ 2 , 

j E dQ 2 

where E is the energy of the parton and Q its virtual mass. The time a parton needs to 
reduce its virtuality from the starting scale Qi ~ E to the hadronization scale Qo — 1 GeV 
can then be estimated as r evo = f dr m — ~. Thus, even though a high energy parton 
with large virtuality will reduce its virtuality rapidly, the overall lifetime is considerable 
and of the order of several fermi. 

Furthermore, the plasma life time r c ~ t (T /T c ) 3 (corresponding to the cool-down from 
the initial temperature to the critical temperature) may be long and of similar magnitude. 
Consequently, the parton shower overlaps in time with the plasma phase. Therefore, split- 
ting and scattering processes have to be treated in a common framework. During the time 
a fast parton reduces its virtuality by building up the parton shower, it experiences scat- 
terings with gluons of thermal mass m s in the plasma. Therefore, we include a scattering 
term S(x, Q 2 ) into the DGLAP evolution for the fragmentation function D m (x, Q 2 ) in the 
medium. The resulting evolution equation combines radiation and scattering: 

dD m {x,Q 2 ) a s (Q 2 ) f 1 dz 



[ ^P(z)D m (*,Q 2 )+S(x,Q 2 ). (5) 

J x z \z / 



d In Q 2 2tt 

For the construction of the scattering term we estimate the relative importance of 
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scatterings with the help of the scattering mean free path A = (rwr) 1 by 

dr _ E± da 2 dtf 

d\ ~ Q2 n dql dQ± Q 2 • {b) 

Two contributions to scattering have to be taken into account: A gain term for scattering 
from a higher energy fraction to the given energy fraction x and a loss term for scattering 
from x to lower energy fractions. Consequently, the energy E in of the incoming parton 
(before the scattering takes place) is different in the gain and the loss term (c.f. Eq. (J7])). 

Scattering makes the fast parton lose energy, which is absorbed as recoil energy q 2 _/(2m s ) 
by the plasma parton. In our treatment, we allow only soft scattering in the t-channel, 
i.e. small relative transverse momentum q\ ~ m|,, such that the change of longitudinal 
momentum fraction Ax = y — x = q\j (2m s E) is small. The Debye mass uid of the gluon 
and the thermal mass m s are related as 2m? s = m 2 D [IE] and mp 3T |34j. In contrast 
to splitting processes, which lead to a multiplicative change of the energy fraction of the 
parton, we model soft scattering processes by shifts to smaller energy fractions. 

Scattering is included in the evolution equation in a similar way as radiation. The 
"scattering probability" Eq. ([6]) is folded with the fragmentation functions and gain and 
loss term are subtracted from each other. The additional scattering term in the evolution 
equation has the following form 

S(x,Q 2 ) = %n [ dw [ dql^(wD m (w,Q 2 )-xD m (x,Q 2 )) 
Q Jx J dq ± 

xsfw-x--^-] . (7) 
We then expand in powers of q 2 _/(2m s E) and drop terms of second order and higher: 



2m s Q 2 J dqY 1 ^ \ v ' ~* ' dx 



2m s Q 2 K D(x,Q 2 ) + x-(x,Q 

Such an expansion is meaningful for small momentum transfer, q 2 _/(2m s E) 1 and x + 
q\l (2m s E) < 1 as we expect them for small momentum transfers q\ ~ m D and large jet 
energies. Remarkably, the scattering term does not depend on the jet energy E in this 
approximation. 
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2 Formalism and Modified Evolution Equations 



In the expression above, n denotes the weighted sum of possible scattering channels. 
Depending on the incident parton, these can be quark-quark, quark-gluon and gluon-gluon 
scatterings: 

- do- 16 3 2ixa s (Q 2 ) 2 j (1 + f) = f incident quark 

n dql n^ [6> {ql + mlf \ (2 + §§) = If incident gluon ' l J 

Interestingly, the quantity appearing in the numerator in Eq. (|E|) is the jet transport 
parameter [5] 

q~na(q 2 ± ), (10) 

which describes the mean acquired transverse momentum per unit length (For a recent 
discussion of q see [T8]). 

For an estimate of the numerical value of q associated with this expression, we use 

, 9 16 . N o 9/^9x (1 + I) = ¥ incident quark . , 

(qi)^ml, q c± -((3)T 3 x 27ral(Ql) x {, I , J, ^ • ., . , ■ (") 



7T- 



+ W = if incident gluon 



The scattering term is most relevant at small virtualities Q ~ Qo and consequently we use 
<2s(Qo) to arrive at an upper boundary for q. Explicitly, Eq. ( ITT]) gives q ~ 0.5 GeV 2 /fm for 
a temperature T = 0.3 GeV for RHIC and q = 5.2 GeV 2 /fm for T = 0.5 GeV corresponding 
to LHC. To fit the experimental data with parton energy loss exclusively, we introduce an 
enhancement factor K in the scattering term. Then the scattering term reads 

s ^= K W{ Dm{x ' Q * )+xdJ ^ M2) )- (i2) 

The value of q = 0.50 GeV 2 /fm (K = 1) has to be viewed in the context of scenario 1 
discussed in the introduction which combines parton energy loss with resonance absorption 
in the fading plasma. Pure parton energy loss in scenario 1 gives a too weak suppression 
factor to explain experimental data from central Au+Au-collisions at RHIC. 

In scenario 2 without final state absorption, we need roughly K = 8 or q ~ 4.0 GeV 2 /fm 
to reproduce the observed nuclear modification factor. If parton energy loss dynamics is 
the same at LHC, this estimate would lead to a value as large as q = 41.9 GeV 2 /fm (for 
K = 8). This increase orginates equally from the higher temperature T = 0.5 GeV and the 
larger scattering cross section of gluons. The values for q for both scenarios are summarized 
in Table [TJ 

In our calculation above, we have introduced K as an artifical enhancement parameter. 
The physical motivation for this soft enhancement factor may come from our estimate of 
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q [GeV 2 /fm] 


T = 0.3 GeV 


T = 0.5 GeV 


Scenario 1 (parton energy loss + absorption), K — 1 
Scenario 2 (large parton energy loss), K — 8 


0.5 
4.0 


5.2 
41.9 



Tab. 1 : Relevant values of the jet transport parameter q in the two scenarios for T = 0.3 GeV 
(incident quark) and T = 0.5 GeV (incident gluon). The effective K -factor allows possible 
enhancements of the energy loss parameter q. 



the scattering cross section. The perturbative running coupling may underestimate the soft 
scattering cross section of the parton in the plasma [32j if one assumes that the temperature 
defines the relevant scale. More precisely, evaluating the running coupling at ttT instead 
of Qo may lead to an enhancement factor of K < 3 for RHIC temperature. A factor of 
K = 8 as in scenario 2 cannot be explained. In contrast, for the higher temperature at 
LHC the scales ttT and Qo are approximately equal which makes an enhancement factor 
K unlikely. 

Up to now, the theoretical formalism has focussed on processes at the parton level. In 
order to compare with observable hadrons, one has to account for hadronization. 

3 Modified Fragmentation Functions D£, DJ in the Quark-Gluon 
Plasma 

For a solution of the modified DGLAP equation ([5]), we evolve a given fragmentation 
function D(x,Ql) to the ultraviolet scale Q max , i-e. we numerically follow a way opposite 
to the shower evolution itself. Furthermore, we assume that the fragmentation function 
at the infrared cut-off scale Qo is unchanged and use standard vacuum fragmentation 
functions as initial condition. We use the AKK parametrization [H] for the fragmentation 
of M-quarks and gluons into pions 

D^(x,Ql) = 0.447 x- L58 (l-x) 101 , ^ ,7r (x, Q 2 ) = 429 x 2m (1 - x) 5 ' 82 (13) 

at the scale Qq = 2 GeV 2 in the evolution equations. The numerical solution of the 
modified DGLAP evolution equation ([5]) (with and without scattering term) is calculated 
with the Runge-Kutta method of 4th order. 

At RHIC the fragmentation of quarks is the main source for pion production. To 
investigate the modification of the ti-quark to pion fragmentation function in the medium 
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3 Modified Fragmentation Functions D%,D* in the Quark-Gluon Plasma 




0.0 0.2 0.4 0.6 0.8 1.0 

X 



Fig. 1: Fragmentation functions D^(x,Q 2 ) in medium and vacuum as a function of energy 
fraction x. The plot compares the results for vacuum evolution ( dashed curve) to medium- 
modified evolution with K = 1 (full drawn thin curve) and K = 8 (full drawn thick curve). For 
the calculation, we have used Q m ax = 

20 GeV and T = 0.3 GeV in Eq. 

at RHIC, we have solved Eq. fl5]) for Qmax = 20 GeV and T = 0.3 GeV. Two different 
values for the enhancement parameter K are used: K = 1 and K = 8. 

The results for the gluon to pion fragmentation function are shown in Fig. [TJ The 
overall behaviour of the medium-modified fragmentation function with Q 2 is similar to 
the vacuum results. The scattering term acts in a similar manner as the splitting term 
and "transports" partons from large to small energy fractions. The suppression of leading 
particles is strongest at large x. Following naive expectations, the suppression of the 
medium fragmentation function becomes stronger for larger K. 

For LHC, gluon fragmentation is more important than quark fragmentation. The re- 
sulting gluon-to-pion fragmentation functions for Q max = 100 GeV and T = 0.5 GeV are 
shown in Fig. [2] for both the vacuum case and for the case of a plasma with K — 1 and 
K = 8. Note that the higher temperature at LHC and the larger color charge of the gluon 
lead to much larger values of q (see Tabled]) although the microscopic dynamics (i.e. energy 
loss of a fast parton) is unchanged. The scattering term drastically reduces the number of 
partons with large x. This effect is much more pronounced at LHC. 

Based on the calculation of modified fragmentation functions, we test the consequences 
of the two discussed scenarios and compute the nuclear modification factor Raa- This 
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Fig. 2: Fragmentation functions Dg(x,Q 2 ) in medium and vacuum as a function of energy 
fraction x. The plot compares the results for vacuum evolution ( dashed curve) to medium- 
modified evolution with K = 1 (full drawn thin curve) and K = 8 (full drawn thick curve). For 
the calculation, we have used Q m ax = 

100 GeV and T = 0.5 GeV in Eq. 



experimental observable can be estimated by folding fragmentation functions with the 
differential cross section for the production of the fast parton. For our comparison, we 
use Pythia 6.4 [15] for a LO calculation of this cross section in pp-collisions at \fs = 200 
GeV. We can fit the Pythia output with the form 



dq 2 ± \ 1 GeV 2 J V 
For quark production at RHIC, the parameters have the following values: 

mb 

A ~ 468.706——, B ~ -3.04, C ~ 9.69. (15) 
GeV^ 

For LHC, the cross section for gluon jet production becomes much larger than for quark 
jets for typical values of p± and consequently our approximation of a gluonic evolution is 
realistic. In the parametrization for -jjp-, we change the parameter values for LHC (with 
y/s = 5.5 TeV) as follows: 

mb 

A ~ 524395——, B ~ -3.27, C ~ 10.75. (16) 
GeV^ 
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4 Effect of Prehadron absorption in Hot Resonance Matter 



Neglecting initial state effects, the nuclear modification factor can be written as a ratio 
of the medium fragmentation function folded with the parton cross section to the vacuum 
fragmentation function folded with the same cross section. 



The results for Raa i n our RHIC scenario for K — 8 (full drawn thick curve) and K — 1 ( 
dashed curve) are shown in Fig. [3j In this plot, it is clearly visible that the natural choice 
K = 1 cannot account for RHIC data for Raa (dashed curve), i.e. in scenario 1 parton 
energy loss alone cannot explain the data. But an artificial tuning of K - and consequently 
of q - to K — 8 matches the observed nuclear modification factor. 

From our perspective, scenario 1 can compete with scenario 2, when an additional 
mechanism for suppression, namely the absorption of resonances in the resonance gas is 
included. This additional physics is quite natural, since the lifetime of the plasma and the 
time for the shower evolution are comparable in magnitude and smaller than the maximum 
path length allowed by the size of the hot medium. Therefore when the plasma fades out, 
confinement is no longer an issue and the created resonance gas will allow prehadrons as 
interaction partners with hadronic cross sections. In the next section, we will discuss in 
detail this scenario of resonance absorption. 

For LHC, we also use Eq. ( TT71) to compute the nuclear modification factor for pions. 
The result is shown in Fig. H]for K = 1 and K = 8. The larger temperature at LHC leads 
to an increase of the suppression relative to RHIC. For the case of K = 8, the computed 
Raa is 10~ 3 — 10 -5 which represents a drastic medium effect. Observing such a large effect 
at LHC would give strong support for a picture in which the suppression is caused by 
parton energy loss with a large transport parameter q. 

The other scenario with K = 1, however, would tell a different story: As discussed in 
the next section, we expect the absorption of resonances to be much less important at LHC 
than compared to RHIC. Consequently, the picture of partonic energy loss becomes more 
realistic and the calculation with K = 1 becomes closer to reality. In this way, the nuclear 
modification factor could increase from RHIC to LHC (i.e. less suppression) although the 
energy loss increases. 



Raa(p±) 



J dzdql^D m (z,Q 2 )6(zq ± -p ± ) 




fdzdqi0rDv(z,Q*)5(zq ± - P± ) 



4 Effect of Prehadron absorption in Hot Resonance Matter 



Various time and length scales are involved in the evolution from the quark-gluon plasma 
to the hadronic phase. Here, we collect numerical values for these time scales to argue that 
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Fig. 3: Two scenarios for the nuclear modification factor Raa of pions at RHIC as a function 
of their transverse momentum p±. The full drawn thin curve corresponds to K = 1 in Eq. \12\i 
while the full drawn thick curve is calculated for a scenario with K = 8. The dashed curve gives 
Raa without absorption. For comparison, experimental data from PHENIX \2U$ are also shown. 



absorption plays an important role at RHIC, but is much less relevant at LHC. 
We consider the following time scales: 

• Equilibration time tq of the plasma, which is estimated to be of the order of 0.5 fm 
at RHIC but much smaller (< 0.2 fm) at LHC [2]. 



Life time of the plasma which we estimate from longitudinal Bjorken expansion 

T r M, | — | (1* 



For an initial temperature T = 0.3 GeV at RHIC and T = 0.5 GeV at LHC, we 
find t c = 3.3 fm at RHIC and r c = 6.1 fm at LHC. 

Average time for evolution of the parton from the hard scale Q max to Qo (where 
prehadrons can be formed): 

(19) 



To- 



Qo Qm 



With E = Q max = 20 GeV for RHIC and E = Q max = 100 GeV for LHC (and 
Qo = V2 GeV), we find r cvo = 2 f m for RHIC and r evo = 10 fm for LHC. Note that 
these time estimate apply to hadrons with average x-fraction. 
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4 Effect of Prehadron absorption in Hot Resonance Matter 
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Fig. 4: Two scenarios for the nuclear modification factor Raa of pions at LHC as a function of 
their transverse momentum p±_ . The dashed curve corresponds to K = 1 in Eq. while the full 
drawn thin curve is calculated for a scenario with K = 8. The dashed curve gives Raa without 
absorption. 



• Time in the hot zone r = R/c: For central Au+Au collisions and central Pb+Pb 
collisions, the size of the plasma is almost identical when estimated in terms of the 
size of the nuclei involved: Ra u — 6.9 fm and i?pb = 7.1 fm. 

The hierarchy between time scales is different at RHIC and LHC. For RHIC, we have the 
following situation: After the evolution of the shower to the virtuality Qq, the plasma has 
become a resonance gas and preconfmed states are formed which can be absorbed in the 
resonance gas. For an estimate of this effect, we use an exponential form for the nuclear 
attenuation, averaged over production points x and path-length (determined by its angle 
<pi of inclination) in a 2-dimensional pancake-like system, 

1 1 f 1 

rab / cl 2 x o d0/— ^e~ nresCTlcsll ^ R ' Xo '^ + e - ™ 105 ^ 05 ' 2 ^^ ' 9 ^) (20) 

In this equation, Z x and li are the paths of the two partons which are flying in opposite 
directions in a medium of size R^n or i?pb respectively. The two respective path lengths 
are determined by |xo ± (li + r c )ei\ = R. For the absorption cross section, we use estimates 
about the resonance gas. These excited hadrons have to deconfme and consequently, their 
size has to be larger than in vacuum. Therefore, the hadronic cross section is estimated to 
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be larger than the irir cross section (and similar to the ttN cross section) because of the 
increased size of the resonances. 

a res ~ 30 mb (21) 

The above cross section concerns hot matter. In cold matter, the medium consists of 
unexcited nucleons. Therefore, the formation of the prehadron plays a more important 
role in the determination of the prehadron-nucleon cross section. Indeed, a cross section 
reduced by a factor of 2/3 compared to its vacuum value fits the suppression rates of 
hadrons in deep-inelastic scattering on nuclei [40J. 

For the density of resonances, we take n res ~ [12]. The duration of this resonance 
phase is considerably longer than the plasma phase, since the entropy density has to be 
dumped into a large volume. 

With the estimate given in Eq. ( 12T1) . one finds values of r a b s — 0.35. The final result 
for the nuclear modification factor, R^a = RAA r abs, is shown in Fig. [3] (full drawn thin 
curve). This scenario of combined suppression mechanisms can principally account for the 
large suppression observed in the nuclear modification factor Raa at RHIC. 

The situation at LHC is, however, qualitatively different. For a 100 GeV jet, the 
development of the parton shower takes about 10 fm and we have r c < r evo . Consequently, 
the parton shower will typically not be surrounded by a deconfined medium already at 
virtuality Qi > Qo. Since r c ~ R, we can estimate an average value for Qi by demanding 

* = !-gf- (22) 



- max 



The finite size of the medium can be included by defining two cut-off scales: At Q\ > Q 
scatterings stop and at Qo perturbative splittings stop. For Q < Q < Q 1 , we use vacuum 
evolution equations. 

A finite medium length reduces the suppression factor slightly. For a medium of length 
R = 7 fm, the nuclear modification factor roughly increases from 0.4 to 0.6 roughly in 
the scenario with K = 1. For the second scenario with K = 8 (large q), we find Raa — 
10~ 2 - 10~ 3 instead of 10~ 3 - 10~ 5 . 

After the evolution to Qo, we may also have resonance absorption at LHC. But since 
r cvo ~10 fm is not much smaller than 2R, absorption cannot play a major role. If we let 
r c take the role of r evo in Eq. ( 1201) . we find a value of r a b s — 0.86 from this equation. 

The final result for scenario 1 (K = 1 including absorption) for LHC is shown in Fig. H] 
(full drawn thin curve). Clearly, absorption is not very important at LHC. 
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5 Modified Parton Distribution in a Jet at small x 



5 Modified Parton Distribution in a Jet at small 



X 



In this section, we study the quark/gluon composition of the jet for RHIC/LHC. At small 
x, the presence of the dense and hot medium modifies the peak position of the Hump- 
Backed Plateau in the ln(l/a;)-distribution. As in the previous section, we consider a 
quark evolution equation in a scenario for RHIC and a gluonic evolution equation for 
LHC. Here, we use the same scattering term as in the previous section. For the evolution 
at small x, a key ingredient is soft gluon coherence. Without soft gluon coherence, the 
parton distribution would diverge at small x. In order to suppress gluon radiation at small 
x, we solve the evolution equation with the scale z 2 Q 2 in the parton distribution [25] 

dD(x,Q 2 ) a s (Q 2 ) f l dz 



dlnQ 2 2tt 



J ^P(z)D (f , z 2 Q 2 ) + S(x, Q 2 ). (23) 



For small x, we can reduce the splitting function to P(z) = 2Ca/z for gluons (LHC) and 
P(z) = 2Cf/z for quarks (RHIC) since these terms represent the dominant contribution. 
In the following, we use the splitting function P(z) = 2Cr/z to account for both cases. 
Such an approach is related to the double-logarithmic approximation (DLA) [29J of parton 
eolution. More terms are taken into account in the framework of MLLA and NMLLA 

[EH EH]. 

The scattering term S(x, Q 2 ) has the same form as in section 2. Its construction and 
relation to the jet transport parameter q is discussed in Eq. flH]). Since the thermal gluon 
mass m s oc T, the scattering term is proportional to T 2 /Q 2 which makes it higher twist, a 
feature shared with the formalisms developed in [331 [22]. Here, we rewrite the scattering 
term as 

2 ^2^ 2 ^ , dD(x,Q 2 " 



S(x, Q 2 ) = ea 2 s (Q 2 )— lD(x, Q 2 ) + x " ^ - 1 (24) 
where 

5AK quark jet 

(25) 

12.2 K gluonjet 

is a dimensionless constant used as an abbreviation for the prefactors. The numerical 
values are calculated with the help of Eq. ( TTTT) . 

The parameter K allows us to study the two scenarios discussed in the previous section: 
Scenario 1 (K = 1) includes small parton energy loss and resonance absorption. Scenario 
2 (K = 8) parametrizes a large transport parameter q tuned up to fit the RHIC data with 
parton energy loss alone. For the calculation of the parton distribution at small x it is 
crucial to keep track of the dependence on a s . We study the most essential modifications 
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and solve the evolution equation in Gaussian approximation [25]. This method reproduces 
experimental data quite well. The dependence of the scattering term on a 2 and the in- 
dependence of the scattering kernel preserves the jet multiplicity, but modifies the position 
of the maximum and the width of the parton distribution in ln(l/x). For further analysis, 
we define the Mellin transform of D(x,Q 2 ) by 

d(J,Q 2 )= [ dxx J ~ 1 D(x,Q 2 ). (26) 



In Mellin space we can simplify the scattering term by partial integration: 

T 2 f 1 / f)D \ T 2 

ea 2 — J dxx J - 1 [D{x ) Q 2 )+x°—{x,Q 2 )^ = ea 2 s — (d(J,Q 2 ) - Jd(J,Q 2 )) (27) 



For constant a s , the evolution equation reads after Mellin transformation 
dd(J, Q 2 ) 



a f 1 T 2 
-i / duP(u)d(j,u 2 Q 2 )-ea 2 — (J - l)d(J,Q 2 ). 

27T Jq Q 



dlnQ 2 

We use the following ansatz for the parton distribution in Mellin space with the anomalous 
dimension 7 

rf(J,g 2 )(x^J (29) 
for which we obtain a consistency equation 7 



7(4 «.) = — T rro (j \ - ec **S( 3 - ^ ( 3 °) 

7r J-l + 2j{J, a s ) Q 2 

At J = 1, the contribution from the scattering term vanishes. Since the first moment with 
J = 1 of the parton distribution is the jet multiplicity, the multiplicity will not change. 
This result is natural for a 2 — >■ 2-scattering term. Modifications coming from the scattering 
term are suppressed since it is a higher-twist cx T 2 /Q 2 contribution. Nevertheless it will 
turn out that the medium significantly modifies the parton distribution D(x,Q 2 ). 
From Eq. ([30]) we find for the anomalous dimension: 



, _ . fT ^fl ea 2 T 2 \ , r ea 2 T 2 \ 2 a s C R 

7 (J |0l ) = -(J-l)( i + Y f ) + ^-l)»( r f ? j +^ (3D 

The small- 2 behaviour of the parton distribution is related to its low Mellin moments, 
therefore we can expand the anomalous dimension around J — 1. 

a s C R , , /l ea 2 T 2 \ 1 HhT , 1 1N2 (I ea 2 T 2 ^ 2 
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5 Modified Parton Distribution in a Jet at small x 



This expansion preserves the most essential features of the parton distribution in vacuum. 
Since the medium does not change the dominant behavior of the splitting term, it also 
works in the medium. For running a s , we use [25] 

d(J, Q 2 ) oc exp (j° J, a s {Q' 2 )) J (33) 

For the calculation of the integral, we use the perturbative 1-loop coupling 

1 11- ^ 

a s (Q 2 ) = r^T^i b = a 3 ' A QCD = 250 MeV. (34) 

6 In t?— 47r 

V QCD / 

The number of flavors is set to five since we study jets with energies well below the 
mass of the top quark. A calculation with a s in one loop approximation is accurate enough 
for the vacuum evolution. Note that d(J, Q 2 ) in Eq. ( 1331) has a Gaussian dependence on 
(J — 1) because the Taylor expansion in Eq. ( 132|) stops at second order: 

d(J,Q 2 ) = (7 exp (a + ai(J - 1) + a 2 (J - l) 2 ) . (35) 

The normalization of d(J, Q 2 ) is not predicted by the ansatz Eq. ( 1331) . We fix the normal- 
ization C ~ 0.024 of the multiplicity n{Q 2 ) = Cexp(a ) with LEP data from e + e" [271128]. 
It is possible to use these "vacuum" data since the scattering term leaves the multiplicity 
unchanged. 

The coefficients in Eq. ([3"5]) result from the integration in Eq. ([3"3"j) and depend on 
both virtuality and temperature. They are given as follows: 



a " = (36) 



6 T 2 



' '' 2 Ei (-(i»,(Q 2 )) -1 ) - [Q 2 -» Ql] 



d2 



2b 2 A 2 
1 

246 



15 6 3 V C^Q 4 v ^ y 26 Y 2C R Q 2 

— T 4 / 2 a s (Q 2 ) 3 / 2 1 a s (Q 2 ) 5/2 



e 2 




15 6 7 /2 A4 fQ 



(37) 



17 





vacuum 


medium 


a 


1 /2C R 

b V ™s 





a% 


1 

Aba , 


-^a 2 


a 2 


./^ 1 a^ 2 


i\J 2C R Q 2 " s 



Tab. 2: Leading terms of the Gaussian coefficients ao, a\ and ai of the parton distribution in a 
jet in Mellin space, i.e. d(J,Q 2 ) ~ exp(ao + ai(J — 1) + a2(J — l) 2 )), which determine multiplicity 
(ao), peak position (a\) and width (ai) of the ln(l/ x)- distribution. The first column gives the 
value in vacuum while the second column shows the correction due to the medium with the lowest 
power in ct s . We use a shorthand notation here: One should take a s — > a s (Q 2 ) in every term 
and subtract the same term with Q 2 — > Qq. 

Here, the notation — [Q 2 — > Qq} means that one has to subtract the previous terms, replac- 
ing Q 2 by Qq. This subtraction originates from the lower limit of integration in Eq. f )33|) . 
Ei(— z) is the exponential integral function, defined by 



Ei(-z) = - 




For a better illustration of our results, we also give the lowest-order expansion of the 
medium modification of the coefficient in a s in Table [2j 

The medium modification is solely contained in the modification of the coefficients 
ai. The coefficients a« have the same physical meaning as in vacuum: ao describes the 
(5 2 -behavior of the jet multiplicity n(Q 2 ) = Cexp(ao), a\ gives the peak position of the 
distribution in ln(l/x) (see below) and a<i is related to the Gaussian width. The first term 
for each coefficient represents the vacuum result while the second terms proportional to 
eT 2 /Q 2 represent the modified evolution in Table El 

In vacuum, the coefficients have inverse powers of a s (Q 2 ). Therefore the vacuum mo- 
ments are determined by the upper virtuality Q 2 with a s {Q 2 ) < a s (Q 2 ) ). The terms 
from the medium-modified evolution have positive powers of a s and change the vacuum 
results. The medium corrections are most important at the infrared scale Qq, where 
a s {Ql) > a s (Q 2 ). Therefore, in the medium the peak position is shifted to larger val- 
ues of ln(l/x) while the width decreases. 
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5 Modified Parton Distribution in a Jet at small x 



The well-known MLLA contribution to the peak position 

1 



a MLLA = 1 f ll^e 27^ 

1 " 2b V 3 3iV c 2 ; V3W c tt 



2\-l/2 



(39) 



cannot be deduced from the Gaussian approximation and consequently is added by hand. 
This contribution is subleading but highly relevant for comparison to experimental data. 
Finally, the parton distribution can be calculated by an inverse Mellin transformation. 



D(x,Q 2 



—— / dJ—j— r exp(a — ai(J — 1) + a 2 (J — 1) 
2ntx J x _ ioo x J ' x 



2irix 
C 

2-kx 



dJexp 



a + iJ \ — ai + In 



a 2 J 



(40) 



Performing the Gaussian integral yields 



xD(x,Q 2 ) 



<Q 2 ) 
2y/ira 2 



exp 



4ao 



(41) 



In Fig. [5] we show a comparison of the logarithmic parton distributions xD(x,Q 2 ) in 
a jet in vacuum and medium. As discussed after Eq. (|25|) . we use the values K = 1 and 
K = 8 and compare to the vacuum results. Two different kinematical cases are studied: 
Q max = 20 GeV and T = 0.3 GeV, suitable for RHIC (top), and Q max = 100 GeV and 
T = 0.5 GeV, relevant for LHC, in the bottom. The shaded area in both figures represents 
estimates for the region where particles cannot be distinguished from the plasma anymore 



At RHIC, the scattering term basically shifts the peak position of the quark distribution 
to larger ln(l/x) (i.e. smaller x) due to the energy loss from gluon scattering with plasma 
gluons. The modification of the parton distribution in scenario 1 with K = 1 is rather 
small. In contrast, the parton distribution in scenario 2 with K = 8 differs enormously 
from the vacuum one. Such a strong modification has not been observed in (preliminary) 
experimental analyses [23]. Note that for jets the existence of the resonance matter is 
unimportant. 

For LHC, particle identification above the soft background from the plasma is possible 
over a much wider range. Scenario 2 with K = 8 leads to a strong modification of the gluon 
distribution: The peak position is shifted to larger values of ln(l/x) and the distribution 
becomes broader. In contrast, scenario 1 with K = 1 again leads only to small changes 
compared to vacuum. 
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Fig. 5: Top: Quark distribution xD(x,Q 2 ) in a jet in Gaussian approximation shown for Q 2 = 
(20 GeV) 2 for evolution in vacuum ( dashed curve) and medium-modified evolution with K = 1 
(full drawn thin curve) and K = 8 (full drawn thick curve) in a plasma with temperature T = 
0.3 GeV as a function of ln(l/x). Bottom: Same for gluons with Q = 100 GeV and T = 0.5 
GeV 
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6 Summary and discussion 



6 Summary and discussion 

We have presented two scenarios for jet evolution in nucleus- nucleus collisions. Scenario 1 
considers weak parton energy loss caused by scattering in the plasma. After the evolution 
to a hadronic scale Qo, preconfinement causes hadronic absorption in the resonance gas 
which is left directly after the plasma decay. In contrast, scenario 2 is solely based on strong 
parton energy loss. In both scenarios, the partonic evolution is computed in the framework 
of a modified DGLAP equation which combines radiation and scattering of the fast parton 
in the medium. The corresponding scattering term is of higher twist, S(x, Q 2 ) oc T 2 /Q 2 . 

In scenario 1, we estimate jet transport parameters of q ~ 0.5 GeV 2 /fm at RHIC and 
q ~ 5.2 GeV 2 /fm at LHC. Both parameters have to be understood as average values since 
we do not consider the change of q with decreasing temperature and density of the plasma. 
Enhancement of the QCD estimate q ~ 0.5 GeV 2 /fm by a i^-factor K = 8 in scenario 2 
can also accomodate our model to the RHIC data for Raa- Such an enhanced value may 
correspond to the hypothesis of a strongly interacting quark-gluon plasma (sQGP). 

We do not favour scenario 2, since there are strong evidences for the importance of 
absorption: (i) Scenario 2 fails to describe production of leading hadrons in deep-inelastic 
scattering [21] ; while models including absorption do well [371 |38j HQ] ; (ii) J/\l/ suppression 
in the hot medium can be described with a rather small value of the transport coefficient 
[42] . which agrees well with the pQCD expectations; (iii) Even with the enhanced transport 
coefficient, scenario 2 does not explain the observed strong suppression of open beauty [39] . 
while the effects of absorption make suppression of light and heavy quarks similar [Jl] ■ 

In contrast, we argue that due to the shorter life times of the shower and the plasma at 
RHIC (compared to the extension of the nuclei), the parton can become preconfined during 
its way through the medium and interact as a prehadron with the decaying resonance 
matter. For an estimate of this effect, we use a standard initial density n Tes ~ T c 3 of 
resonances and an interaction cross section of a res = 30 mb. With these values, we obtain 
an acceptable fit to the RHIC data with K = 1. 

The new energy domain of LHC allows to test these two alternatives in two different 
ways. First, we think that due to the similarity of a s (nT) and a s (Qo) for T = 500 MeV, 
K — 1 is also reasonable for LHC. Then, we predict in scenario 1 a slightly larger Raa due 
to a stronger partonic energy loss while hadronic absorption plays a minor role. At LHC, 
the length of the perturbative shower is much longer, leaving almost no time for resonance 
absorption in the medium after preconfinement. Consequently, prehadron absorption does 
not increase the suppression factor significantly at LHC. In contrast to scenario 1, scenario 2 
would lead to a very large suppression, i.e. Raa < 10~ 2 at LHC. Such a strong suppression 
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may be detected within the first few months of the heavy-ion program at LHC. 

Secondly, one can study the distribution of low-x partons in the jet, since at LHC 
particles up to ln(l/x) ~ 5 are measurable above the plasma background. The shape 
of the hump-backed plateau of the medium-modified jet will clearly show whether any 
enhancement of q in our calculation is necessary. 

For our considerations, we have not made use of further experimental data [13J HU H5] 
for several reasons. First, the results from SPS in Ref. [13] and RHIC in Ref. [H] are 
not fully consistent with each other. The absence of suppression at small energies in jB] 
is unexpected in our absorption scenario. However, Ref. jl3] finds a suppression at a 
similar energy in the most central collisions. Furthermore, we mention as an aside that the 
suppression is estimated to set in after a propagation length of 2 fm in [15] which may be 
related to our estimate for the shower time. 

As a final remark, we note that the jet composition in the medium from our calculation 
can easily be distinguished from approaches with medium- induced gluon radiation [46]. 
The enhancement of the splitting functions in these calculations always leads to a much 
higher peak and consequently to a much larger multiplicity |47J. 
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